Low energy atmospheric muon neutrinos in MACRO 



The MACRO Collaboration: 

M. Ambrosio^^, R. Antolini'^, G. Auriemma"'" , D. Bakari^'i^, A. Baldini^^, G. C. Barbarino^^, B. C. Barish'', 
G. Battistoni^'^ R. Bellotti\ C. Bemporad", P. Bernardinii°, H. Bilokon^ V. Bisil^ C. Bloise^ C. Bower*^, 
M. Brigida\ S. Bussinol^ F. Cafagna\ M. Calicchio\ D. Campana^^ M. Carboni^, S. Cecchini^'^, F. Cei^^^^^ 
V. Chiarella^ B. C. Choudhary*, S. Coutu"'', G. De Cataldo\ H. Dekhissi^-", C. De Marzo^, I. De Mitrii'', 
J. Dcrkaoui^'^'', M. De Vincenzi^^, A. Di Credico'', O. Erriquez^, C. Favuzzi^, C. Forti^, P. Fusco^, G. Giacomelli^, 
G. Gianninii^'^ N. Giglietto\ M. Giorgini^, M. Grassi^^, L. Gray"^, A. Grillo^, F. Guarino^^, C. Gustavino'^, 
A. Habig'^, K. Hanson-'^^, R. Heinz^, E. larocci^'-^, E. Katsavounidis^, I. Katsavounidis*, E. Kearns'^, H. Kim^, 
S. Kyriazopoulou^, E. Lamanna"'^ C. Lane^, D. S. Levin^^, P. Lipari^'^, N. P. Longley'*'*, M. J. Longo^^, 
F. Loparco^ F. Maaroufi^'i"^, G. Mancarella^o, G. Mandrioli^, A. Margiotta^, A. Marini^, D. Martello^", 
A. Marzari-Chiesai*', M. N. Mazziotta\ D. G. Michael"*, S. Mikheyev^^'^'f, L. Miller®'", P. Monacelli^, T. Montaruli^ 

M. Montenois, S. Mufson*, J. Musser^ D. Nicolo".'', R. Nolty*, C. Orth^, G. Osteriai2, M. Ouchrif^'i^^ 
O. Palamara^ V. Patera^'-'', L. Patrizii^, R. Pazzi^^, C. W. Peck"*, L. Perrone^o, S. Petrera^ P. Pistilli^®, V. Popa^''' 
A. Raino\ J. Reynoldson"^ , F. Ronga*^, C. Satriano^''''", L. Satta'^'-'', E. Scapparone"^, K. Scholberg^, A. Sciubba'^'-'', 
P. Serra^, M. Sioli^, G. Sirri^.^, M. Sitta^^, P. Spinellii, M. Spinetti^ M. Spurio^'*, R. Steinberg^, J. L. Stone^, 
L. R. Sulak^, A. Surdo^o, G. Tarle V. Togo^, M. VakiRis, E. Vilela^, C. W. Walter^'* and R. Webb^^ 
1. Dipartimento di Fisica dell'Universita di Bari and INFN, 70126 Bari, Italy 
2. Dipartimento di Fisica dell'Universita di Bologna and INFN, 40126 Bologna, Italy 
3. Physics Department, Boston University, Boston, MA 02215, USA 
4- California Institute of Technology, Pasadena, CA 91125, USA 

5. Department of Physics, Drexel University, Philadelphia, PA 19104, USA 

6. Laboratori Nazionali di Frascati dell'INFN, 00044 Frascati (Roma), Italy 

7. Laboratori Nazionali del Gran Sasso dell'INFN, 67010 Assergi (L'Aquila), Italy 
8. Depts. of Physics and of Astronomy, Indiana University, Bloomington, IN 47405, USA 
9. Dipartimento di Fisica dell'Universita dell'Aquila and INFN, 67100 L'Aquila, Italy 

10. Dipartimento di Fisica dell'Universita di Lecce and INFN, 73100 Lecce, Italy 

11. Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA 
12. Dipartimento di Fisica dell'Universita di Napoli and INFN, 80125 Napoli, Italy 

13. Dipartimento di Fisica dell'Universita di Pisa and INFN, 56010 Pisa, Italy 
14. Dipartimento di Fisica dell'Universita di Roma "La Sapienza" and INFN, 00185 Roma, Italy 

15. Physics Department, Texas A&M University, College Station, TX 77843, USA 
16. Dipartimento di Fisica Sperimentale dell'Universita di Torino and INFN, 10125 Torino, Italy 
17. L.P.T.P., Faculty of Sciences, University Mohamed I, B.P. 524 Oujda, Morocco 
18. Dipartimento di Fisica dell'Universita di Roma Tre and INFN Sezione Roma Tre, OOI46 Roma, Italy 

a Also Univorsita della Basilicata, 85100 Potenza, Italy - b Also INFN Milano, 20133 Milano, Italy 
c Also Istituto TESRE/CNR, 40129 Bologna, Italy - d Also Scuola Normale Superiore di Pisa, 56010 Pisa, Italy 
e Also Universita di Trieste and INFN, 34100 Trieste, Italy - / Also Dip. di Energetica, Universita di Roma, 00185 Roma, Italy 
g Also Institute for Nuclear Research, Russian Academy of Science, 117312 Moscow, Russia 
h Also Institute for Space Sciences, 76900 Bucharest, Romania - i The Colorado College, Colorado Springs, CO 80903, USA 
I Also Dept. of Physics, Pennsylvania State University, University Park, PA 16801, USA 
m Also Dipartimento di Fisica dell'Universita della Calabria, Rende (Cs), Italy 
n Also Dept. of Phys. James Madison University, Harrisonburg, VA 22807, USA 
* Corresponding author, e-mail: SPURIO@BO.INFN.IT 

We present the measurement of two event samples induced by atmospheric z/^ of average energy 
Ev ~ 4 GeV . In the first sample, the neutrino interacts inside the MACRO detector producing 
an upward-going muon leaving the apparatus. The ratio of the number of observed to expected 
events is 0.57 ± 0.05stat ± Q.QQsyst ± 0.14tj,eor with an angular distribution similar to that expected 
from the Bartol atmospheric neutrino flux. The second is a mixed sample of internally produced 
downward-going muons and externally produced upward-going muons stopping inside the detector. 
These two subsamples are selected by topological criteria; the lack of timing information makes it 
impossible to distinguish stopping from downgoing muons. The ratio of the number of observed to 
expected events is 0.71 ± O.OSstat ± 0.07syst ± O.lStheor • Using the ratio of the two subsamples (for 
which most theoretical uncertainties cancel) we can test the pathlength dependence of the oscillation 
hypothesis. The probability of agreement with the no-oscillation hypothesis is 5% . 

The deviations of our observations from the expectations has a preferred interpretation in terms 
of oscillations with maximal mixing and Am^ ~ 10~^ 10~^ eV^ . These parameters are in 
agreement with our results from upward throughgoing muons, induced by i/^j of much higher energies. 
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The results from several underground detectors which 
measure the flux of atmospheric neutrinos give strong in- 
dication that v^^s oscillate into neutrinos of another type 
10-^. Fully-contained and partially-contained neutrino- 
induced events observed in underground detectors come 
from neutrinos of energy ~ 1 GeV. The flux of atmo- 
spheric neutrinos of several tens of GeV can be inferred 
from the measurement of neutrino-induced upward-going 
muons that traverse the entire detector {up-throughgoing 
muons). The hypothesis of neutrino oscillations, with 
best-fit parameters sin^ '^-O^ix ^ 1 a-nd Am^ in the range 
of a few times 10"'^ eV^, can explain the observed anoma- 
lies both in the ratio of contained Vfj, to Uf, events (Super- 
K, Soudan 2) and in the zenith angle distribution of up- 
throughgoing muons (MACRO, Super-K). 

The MACRO detector measures both the high energy 
(median energy 50 GeV) and few GeV energy at- 
mospheric neutrino fluxes. In Ref. the interpreta- 
tion of the data in terms of v^^ oscillations came from a 
deficit and from an anomalous zenith angle distribution 
of the observed up-throughgoing muons originating from 
interactions in the rock below the detector. Here we 
report on the measurement of the flux of lower energy 
{Ey ~ 4 GeV) atmospheric neutrinos through the detec- 
tion of interactions inside the apparatus (yielding par- 
tially contained upgoing and downgoing muons) and by 
the detection of externally produced upward-going muons 
stopping inside the detector 

MACRO 1^ is a large area, modular tracking detector 
located in Hall B of the Gran Sasso Underground Labo- 
ratory in Italy, with an average rock overburden of 3700 
hg/cm^. It is a rectangular box, 76.6 m x 12 m x 9.3 m, 
divided longitudinally into six supermodules and verti- 
cally into a lower part (4.8 m high) and an upper part 
(4.5 m high). The active detection elements are planes of 
streamer tubes for tracking and liquid scintillation coun- 
ters for fast timing. The lower half of the detector is 
filled with streamer tube planes alternating with trays of 
crushed rock, which provide most of the 5.3 kton target 
mass for partially-contained neutrino interactions. The 
upper part is hollow and contains the electronics racks 
and work areas. There are 10 horizontal streamer tube 
planes in the bottom half of the detector, and 4 planes 
on the top, all with wire and 27° stereo strip readouts. 
Six vertical planes of streamer tubes cover each side of 
the detector. The intrinsic angular resolution for muons 
is between 0.1° and 1.0° depending on the track length. 
The scintillator system consists of three widely-separated 
layers of horizontal boxes, and on each vertical side of 
the detector a layer of vertical boxes inserted between 
the streamer tubes. The time (position) resolution for 
muons in a scintillator box is about 500 ps (~ 11 cm). 
The direction of the muons passing across MACRO is 
determined by the time-of-flight between two layers of 
scintillation counters. 

The results presented in this letter come from 4.1 live 
years of data taking with the full detector, from April 
1994 to February 1999. 
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FIG. 1. Event topologies induced by neutrino-'interactions 
in or around MACRO. lU = Internal Upgoing /i; ID= Internal 
Downgoing /x; UGS= Upgoing Stopping fi; Up throughgoing 
— upward throughgoing fi. The the black circles indicate 
the streamer tube hits, and the black boxes the scintillator 
hits. The T.o.F. of the muon is measured for the lU and Up 
throughgoing events. 

About 33 X 10^ downgoing muons were collected, and 
were used to monitor the detector efficiency, the running 
conditions and the acceptance. The trigger rate due to 
downgoing muons is 0.3 Hz. The trigger efficiency for 
each scintillation counter and for the streamer tubes was 
monitored over the data taking period using the down- 
going muons. 

Two samples of atmospheric muon neutrinos in the 
few-Gel^ energy range are measured. In the first sam- 
ple (up partially-contained or /[/=Internal Upgoing /i 
events) there are (mainly) events induced by charged cur- 
rent (CC) interactions of upgoing v^^ inside the lower part 
of MACRO. An upgoing muon is produced, which crosses 
two scintillation layers (Fig. |l|), so that the measurement 
of the direction is made through time-of-flight. 

The second sample is a mix of upgoing and downgoing 
events. The partially contained downgoing events (down 
partially-contained or /Z?=Internal Downgoing /x) are in- 
duced by downgoing v^, interacting in the lower part of 
MACRO. The upgoing stopping events ( UGS = Upward 
Going Stopping muons) are induced by interactions of 
upgoing below the detector yielding an upgoing muon 
which stops inside the detector. Both the down partially- 
contained and the upgoing stop events cross only the bot- 
tom layer of liquid scintillation counters (see Fig. |l|) and 
are identified by means of topological criteria. The lack 
of timing information makes it impossible to distinguish 
between the two subsamples. Fig. ^ shows the parent 
neutrino energy distribution from a Monte Carlo calcula- 
tion for the three event topologies detectable in MACRO. 
The energy spectrum and the median energy of the two 
samples presented in this letter are almost the same. 
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FIG. 2. Monte Carlo simulated distribution of the parent 
neutrino energy giving rise to the three different topologies of 
events detectable by MACRO. The distributions are normal- 
ized to one year of data taking; the analysis cuts are included. 



The identification of lU events is based on topological 
criteria and time-of-flight measurements. The main re- 
quirement is the presence of at least two hit scintillator 
clusters, respectively in the center layer and in the upper 
part of the apparatus (see Fig. |l]). This is the expected 
topology for a neutrino interacting in the lower detec- 
tor and producing an upward-going muon with enough 
energy to exit the apparatus. It is also the topology of 
the much more numerous downgoing muons stopping in 
the lower detector. Scintillation timing allows the sepa- 
ration of the two classes of events. Moreover, the scin- 
tillators are required to match a streamer tube track re- 
constructed in space by our standard track-finding al- 
gorithms p. For lU candidates, the lowest point of a 
track (the starting point) must be inside the apparatus 
as a condition for the containment of the i/^ interaction 
vertex. To reject fake partially-contained events entering 
from a detector insensitive zone, the extrapolation of the 
track in the lower part of the detector must geometrically 
cross and not fire at least one scintillator layer and one 
streamer tube plane, or at least three planes of streamer 
tubes. These conditions were tuned on Monte Carlo sim- 
ulated events, including evaluation of detector inefficien- 
cies. Other cuts are applied to reject background events 
from downgoing atmospheric muons. They are related 
to the goodness of the geometrical agreement between 
scintillator hits and the streamer track, to the proper op- 
eration of the scintillation counters and to the quality 
of the time measurement. The measured 1//3 distribu- 
tion after all analysis cuts (including the requirements of 
vertex containment) is shown in Fig. ^. The measured 
muon velocity /3c is evaluated with the sign convention 
that upgoing muons have 1/f] ^ —I. A total of 121 events 
survive in the range —1.3 < 1//3 < —0.7, which is taken 
as the range of lU signal. 

From the time distribution of Fig. one expects some 



background events; they are mostly due to wrong time 
measurements or secondary particle hits, yielding an al- 
most flat l/f3 distribution. The fit of the distribution 
in the range —6.0 < 1//3 < —0.3 to a gaussian plus a 
straight line gives an estimated background of 5 events 
in the signal region. After background subtraction, we 
have 116 up partially-contained events. 
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FIC. 3. The 1//3 distribution of partially-contained events 
after all software cuts. 121 events are in the lU event signal 
region, —1.3 < 1//3 < —0.7, with 5 estimated background 
events. 

The identification of ID + UGS events is based on 
topological criteria. The candidates have a track start- 
ing (ending) in the lower apparatus, and crossing the 
bottom detector face. The track must also be located 
or oriented in such a way that it could not have entered 
or exited undetected through insensitive zones in the ap- 
paratus. Events with scintillator hits outside the bot- 
tom layer, or with the reconstructed track pointing to a 
detector insensitive zone between modules, are rejected. 
The standard muon tracking procedure B is based on 
at least four aligned hits in the streamer tubes. This 
corresponds to a minimum traversed detector thickness 
tv ~ 200 g/cm? (standard sample). A dedicated retrack- 
ing procedure was applied to all remaining events. The 
retracking procedure requires at least three streamer tube 
hits {t^ ~ 100 g/cm^), ahgned with respect to a fired 
scintillation counter in the bottom layer. The number of 
retracked events is ~ 5 x 10^ (to be compared with the 
~ 32 X 10^ standard muon tracks); these events follow 
the same analysis as the standard sample. 

After the software cuts, 879 events survive. Some of 
them are tracked incorrectly (mostly by the retracking), 
or are bending downgoing muons, entering from a detec- 
tor insensitive zone. Due to the bending, only a fraction 
of the streamer tube hits are used by the tracking algo- 
rithms. In order to reject these fake candidates, we made 
a double scan with the MACRO Event Display [§. To 
eliminate any bias from the scan procedure, and to evalu- 
ate the absolute and relative scanning efficiencies, Monte 
Carlo (MC) simulated events (described below) were ran- 
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domly injected into the data sample before the scanning 
stage. Two physicists independently scanned the merged 
sample. At the end of the scan, 200 events in the real 
data (106 ID + UGS candidates are retracked events) 
are accepted as upgoing stopping or partially contained 
downgoing muons. 97% of the real events selected by one 
physicist were also selected by the other. 

Downgoing muons which pass near or through the de- 
tector may produce low-energy, upgoing particles, which 
could simulate neutrino-induced upgoing muons if the 
downgoing muon misses the detector. This background 
has been evaluated using a full simulation, based on our 
measurements Q]. The background is 7 ± 2 events. A 
second background source could arise from atmospheric 
muons and detector inefficiency. Using a simulated sam- 
ple of 10^ atmospheric muons, which includes measured 
detector inefficiencies, no events were selected by the 
above described procedure. After background subtrac- 
tion, 193 events represent the down partially-contained 
plus upgoing stopping signal. 

The expected number of neutrino-induced events was 
estimated from GMACRO §, a GEANT-based @ full 
MC detector simulation. The Ve and interaction rates 
have been computed using the atmospheric neutrino flux 
of the Bartol group ||] and the neutrino cross sections of 
Ref. 0. In this cross-section model, the contributions 
of the exclusive channels of lowest multiplicity (quasi- 
elastic and single pion production) are calculated sep- 
arately from deep inelastic scattering (DIS). The DIS 
contribution to the vN cross section was computed us- 
ing the GRV-LO-94 |jll| parton distribution functions. 
Using these neutrino fluxes and cross sections for E^, > 
300 MeV, we expect a total (CC + NC) interaction rate 
of 71.5 /kton-y [ve+Ve) and 148.1 /kton-y {v^+v^). Two 
simulated samples have been generated, because of the 
different vertex locations for the lU and the ID + U GS 
events; the simulated events were processed with the 
same analysis chain as the data. 

For the lU events, a sample of 10^ interactions inside 
the apparatus was generated (equivalent to 85.9 years 
live time). The simulation indicates (see Table 1) that 
87% of detected lU events come from j/^ charged cur- 
rent (CC) interactions, 9% from ly^ CC and the remain- 
ing fraction from neutral current (NC) interactions. Due 
to detector inefficiencies and analysis algorithm failures, 
some neutrino- induced events originating in the rock sur- 
rounding the detector are expected to contribute to the 
selected sample of up partially-contained events (upward- 
throughgoing /i's appearing as partially-contained). The 
vertex containment requirements reduce this background 
to about 1%, evaluated using a simulated sample of up 
throughgoing muons |3|. The fully-automated selection 
gives a total number of 202 expected events in the lU 
event signal region, —1.3 <!/ (3 < —0.7 for 4.1 y of live- 
time . 

For the ID -\- UGS events, 1.16 x 10^ neutrino interac- 
tions were simulated in a larger volume (including the ex- 
perimental hall, the detector and the surrounding rock). 



The generated events correspond to a live time of 31.1 y. 
The dimensions of the interaction volume (with 13 m of 
rock below the detector, and a total rock mass of 165 kton 
plus 5.3 kton of the apparatus itself) were chosen to re- 
duce to less than 1% the number of i^^-induced stopping 
muons produced outside that volume. The 2199 events 
which survived the software selection for the ID + UGS 
were merged with the real events which passed the same 
software selection, and visually scanned. After the scan 
procedure 2074 (=94.3%) of the simulated and recon- 
structed events were accepted, together with the 200 real 
events. The expected rate is 273 events in 4.1 y live- 
time. In Table | we give the main features of the lU and 
ID + UGS simulated samples (expected rate, percentage 
of CC interactions, median parent neutrino energy, 
and fraction of events induced by upgoing neutrinos). 

The number of detected lU events in the 4.1 y live- 
time is 116, while the expected number is 202. The ra- 
tio of the measured to the expected events is Riu = 
{^)iu = 0.57 ± 0.05stat ± OMsyst ± O.Utheor. For 
the ID + UGS sample, 193 events are detected while 273 
are expected. The ratio is Rid+ugs = 0.71 ± 0.05stat ± 
0.07 sy St i O.lStheor- Each data set is (within errors) con- 
sistent with a constant deficit (43% for the lU sample, 
29% for the ID-f-UGS) in all zenith bins compared to 
the Monte Carlo expectations assuming no oscillations. 
Fig. ^ shows the zenith angle distributions for the two 
measured data sets and for the Monte Carlo simulations. 
Due to the analisys cuts and to the apparatus acceptance, 
there are no events in the last bin of the two distributions. 
The expectations are affected by a systematic theoretical 
error due to the uncertainties regarding the atmospheric 
neutrino flux and the neutrino cross sections. At present 
there is no unique and reliable estimate of the total theo- 
retical uncertainty for the rate calculations. Each exper- 
imental group, and for each event category, has its own 
way to estimate it. For this analysis we conservatively 
estimate 20% for the flux and 15% for the cross section, 
which add in quadrature to an error of 25%. For our 
high energy events [|j we quoted 17%, while in the re- 
cent SuperKamiokande analysis of neutrino-induced 
stopping muons 22% was quoted . 

Our measured value of Riu = 0.57 is quite far from its 
expected value of unity. If we ignored theoretical errors 
(i.e. if we assumed the flux and cross section as we mod- 
eled them were accurate), an experiment with our statis- 
tical and experimental uncertainties would only fluctuate 
so far from unity with probability 2.5 x 10~^. However, if 
we take the 25% theoretical error into account, the prob- 
abihty becomes 6.5%. Rid+ugs — 0.71 also differs from 
unity, though not as significantly as Riu- 

If the observed deficit were due only to an overall the- 
oretical overestimate of the neutrino flux or cross sec- 
tions, one would expect Rju ~ Rid+ugs (small differ- 
ences would remain due to residual geomagnetic effects) . 
Furthermore, the theoretical uncertainties largely cancel 
if the ratio lU/ ID + UGS between the measured num- 
ber of events is compared with the expectation. The 
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partial error cancellation arises from the almost identi- 
cal energy spectra of parent neutrinos for the two sam- 
ples of events; we evaluated the remaining error as 5%. 
The experimental systematic uncertainty for the ratio 
is estimated at 6%. The measured ratio is 



lU 



ID+UGS ~ 

0.60 ±0.07 Stat, while the expectation without oscillations 
is 0.74 ± O.Oisys ± O.OAtheo- The probability to obtain a 
ratio so far from the expected one is 5%, almost inde- 
pendent of the neutrino flux and neutrino cross sections 
used for the predictions. 
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FIG. 4. Measured distributions in the cosine of the zenith 
angle O for the a) ID + UGS and h)IU events (black points 
with error bars). The solid lines are the Monte Carlo predic- 
tions assuming no oscillations. The dashed lines are the ex- 
pectations for i^fi Ut oscillations with 
and maximal mixing. 

We investigated if the observed discrepancies between 
data and expectations could be explained by possible sys- 
tematic effects. The detector mass is known to ±5%. 
The uncertainty for the detector acceptance was esti- 
mated by comparing the shape of the zenith distribution 
of downward-going muons stopping inside the detector 
with a MC expectation based on the known rock over- 
burden: the two distributions agree within 6%. Other 
uncertainties arise from the live-time estimate (3%), the 
effective containment of the interaction vertex depend- 
ing on the simulation of the detector response to internal 
neutrino interactions (4%) and the background subtrac- 
tion (4%). Adding all these contributions in quadrature 
yields our quoted experimental systematic uncertainty of 
10%, too small to account for the observed discrepancy. 
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FIG. 5. Distribution of the amount of material {t in 
g cm~^) in the detector traversed by lU and ID + UGS 
events. The data points (black points with error bars) are 
compared with the Monte Carlo prediction and no oscillations 
(solid histogram). The expectation for oscillations 
with maximal mixing and Am^ = 2.5 X 10 ^eV^ is shown by 
the dashed histogram. 

The number of expected events was also evaluated us- 
ing the NEUGEN neutrino event generator ||l3| (devel- 
oped by the Soudan and MINOS collaborations) as input 
to our MC simulations. The NEUGEN generator pre- 
dicts ~ 6% (5%) fewer lU (ID-fUGS) events than our 
default generator pO[ | , well within the estimated system- 
atic uncertainty for neutrino cross sections (^ 15%). 

Our data disfavor the no-oscillations hypothesis re- 
gardless of overall normalization; they are consistent with 
neutrino oscillations (i/^ disappearance) with maximal 
mixing and Am^ ^ (1^20) x 10~^ eV^. As a "test point" , 
we use the best-fit parameters from our high-energy anal- 
yses Am2 = 2.5 X 10^3 eV^ and sin^2B^i^ = 1. 
The predicted numbers of events and the angular distri- 
butions are indicated by the dashed histograms in Fig. 
0; they are in good agreement with the measured data. 



For Ar 



(1 20) X 10 eV , upgoing neutrinos 



(which induce lU and UGS events), which travel thou- 
sands of kilometers through the Earth, are reduced by 
50%. Almost no reduction is expected for downgoing 
partially-contained muons. In this scenario, and for a 
pure CC interaction sample, the expected event rate 
is 1/2 of the lU and 3/4 of the ID + UGS predictions 
without oscillations. The predicted reduction for upgo- 
ing is less than 1/2 because of the I'e and NC event 
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contaminations. Our data disfavor Am^ > 10~^ eV'^, for 
which the ID events are also reduced; both the ID+UGS 
and lU event rates are ^ 1/2 of the no-oscillations expec- 
tation. We also disfavor Am^ < 10"'^ eV"^, for which the 
shape of the angular distributions (Fig. ^) is modified. 

Assuming oscillations (with the "test point" parame- 
ters) 115 up partially-contained and 202 down partially- 
contained plus upgoing stopping events are expected. For 
the lU events, the reduction from the no oscillations hy- 
pothesis is 0.57, to be compared with the measured value 
of Riu 0.57. For the ID + UGS events, it is 0.76, to be 
compared with Rjd+jjcs = 0.71. The quoted numbers 
use our default normalization. Recent flux calculations 
psf suggest that the Bartol flux which we use may be too 
high (though within the quoted theoretical error). As far 
as the event rates are concerned, a lower normalization of 
the flux can still be partially compensated at low energies 
by different interaction cross sections for neutrinos. 

The event distributions as a function of the muon path- 
length inside the detector have also been studied, as an 
independent consistency check. In Fig. || the data (black 
points) are compared with the MC expectation (solid 
lines for no neutrino oscillations; dashed lines for oscilla- 
tions at our test point). The shapes are similar, and the 
data prefer the reduced normalization of the oscillation 
prediction. 
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FIG. 6. Allowed contours at 90% C.L. for i/f^ — > Ut oscil- 
lations obtained by combining the low energy neutrino events 
{lU and ID + UGS) using the prescription of @. The 90% 
C.L. contour and the best fit point from our high-energy anal- 
ysis is also shown as a dashed line. 

We estimated the most likely values of Am^ and 
sin'^26mix using a method for data and Monte Carlo 
for the data of Fig. ^. The was computed with ten de- 
grees of freedom: the histograms (2x4 bins, normalized 
so as to contain only distribution shape information) , the 
ID+UGS '"^tio and the overall normalization. The sta- 
tistical and systematic errors are added in quadrature; 



the systematic uncertainty is 10% in each bin of the an- 
gular distributions, 5% for the ratio, and 25% for the 
normalization. Fig. ^ shows the 90% confidence level 
region, based on the application of the MC prescriptions 
of Ref . [|l6| on a {si'n?26mix, Arn^) grid. The expected 
flux $^3^ for a given point of {sin'^26„iix, Arri^) in the 
grid is obtained by weighting each simulated event with 
its surviving probability P{Vf_i — > v^j,) = 1 — sin^29mix ■ 
sm^(1.27 • Am^L/E^) in that bin. The maximum of 
the probability (97%) occurs at sin'^29mix = 1-0; 
this probability is almost constant in the interval 
Am^ = (1 20) X 10"'^ eV'^. In the region of the maxi- 
mum, the flux normalization factor is 1.02 in both data 
sets ( i.e. the data are 2% higher than the oscillated 
predictions with our normalization). 

In conclusion, we presented measurements of two sam- 
ples of events induced by relatively low-energy neutrinos 
(El, ~ 4 GeV) interacting in MACRO or in the sur- 
rounding rock. The neutrinos originate from cosmic ray 
interactions in the upper atmosphere above the detec- 
tor (downgoing events) or up to ~ 13000 km away on 
the opposite side of the Earth (upgoing events). The 
ratio of the number of observed to expected events (no 
oscillations) is 0.57 ± 0.05stat ± Q-O&syst ± ^-^^theor for 
the lU events and 0.71 ± 0.05siat ± Q-Olsys ± G-l^theo for 
the ID -I- UGS. Within statistics, the observed deficits 
are uniform over the zenith angle. From the ratio of lU 
to ID + UGS, the probability that there is an overall 
reduction in the number of neutrino-induced muons is 
5%. The hypothesis of muon neutrino oscillations ex- 
plains the different deficits in lU and ID + U GS events 
with higher probability. The large theoretical uncertain- 
ties for the neutrino fiux and cross sections is dominant in 
each data set; the ratio of the two low energy samples is 
dominated by statistical uncertainties. The regions with 
Am^ > 3 X 10^"* eV^ and siri?'29mix > 0.25 are allowed at 
90% C.L. The best region corresponds to maximal mixing 
and Am^ = (l-=-20) x 10~^ eV'^. This result confirms the 
scenario proposed by the measurement of higher-energy 
neutrino- induced muons by MACRO as well as by 

other experiments 0J^, all of which favor the oscilla- 
tion hypothesis with maximal mixing and Arn^ of a few 
times 10~^ eV'^. 





Rate 


/^M CC 








(ev/y) 


(%) 


(GcV) 


(%) 


lU 


49.3 


87 


4.2 


94 


ID+UGS 


66.7 


87 


3.5 


51 



TABLE I. Summary of the properties of simulated lU and 
ID + UGS events. From column 2 to 5: expected event rate 
(no neutrino oscillations); percentages of events induced by 
CC interactions; median parent neutrino energies (5% of 
events with the highest /lowest energies were excluded); per- 
centages of events induced by upgoing neutrinos. 
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